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I. INTRODUCTSON 


When optical radiation passes through the atmosphere the 
amplitude and phase characteristics of the beam are altered 
by the. inhomogeneity of the index of the refraction along 
the propagation path. These effects limit the performance 
of laser communication and radar systems. It is, therefore, 
important to have a theory for predicting the nature of the 
propagation of the light beam in a random medium. 

Dr. Ze'evi [Ref.1] in his doctoral thesis developed a 
weighting function for the scintillation of a spherical wave 
over a folded path. To verify Ze'evi's function many exper- 
iments were attempted by Speer and Parker [Ref.2], 
Constantine [Ref.3], Flenniken [Ref.4], and Henry [Ref.5]. 

Essentially our work was divided into three parts. The 
first part was the characterization of turbulence in the 
turbulence chamber and the second part was the measurement 
of the experimental weighting function for the folded path. 
The third part was to compare the experimental results to 
the Ze'evi's weighting function and to describe the results 
of the effects of the geometric optical model of turbulence. 

In this paper we concentrated on the third part. The 
first and the second parts will be discussed in detail in 
the papers of Mr. Lee [Ref.6] and of Mr. Kim [Ref.7] 
respectively. 

We first review the development of Dr. Ze'evi's 
weighting function for the folded path and discuss briefly 
the experimental procedure. Then, we discuss the effects of 


the geometrical-optic model of a turbulent medium. 
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Pie lhaeOnn tT TCAL BACKGROUND 


Here, we will first outline the theory of the index of 
refraction fluctuations and then review the log-amplitude 


variance and finally Ze'‘evi's folded path theory. 


ma. INDEX OF REBRACTION 


In optical propagation through the atmosphere, the fluc- 
tuations of refraction index are sufficient to produce beam 
wander, beam spread, ANGI neenstty flmetuations . The 
refractive index is assumed to be a real function of the 
temperature and the pressure, and has the following 
form [Ref.8]: 





ne Oe 4, e 
nin) =10+ Spur) (1+ ar) * io (eqn 2.1) 
where n(r) = real refraction index at a point r, 
T(r) = Temperature ( Kelvin ), 
P(r) = Pressure Guia liabar 


aN 


Beam wavelength ( micrometer ). 


The dependence of the refractive index of the atmosphere on 
the change in P andT is given by the differential of 


equation 2.1 





| mie de ale Bie Taeae 
= Pee ae a ‘ ie eqn 2.2 
dn= 776 = ( D =e xD dx ite (eq ) 

Since the pressure fluctuations are relatively small, 
the fluctuations of the refraction index are governed by 
temperature fluctuations. An accepted measure of the 


strength of turbulence in the atmosphere is the refractive 


Ll 


index structure parameter, C,*, which is proporttonal spose 
fluctuations of the refractive index. 
The structure function of the refractive index “ine 


Kolmogorov inertial subrange is: 


D,(ny = <[ni ) > 1G (eqn =e 
where Il1<<| r - fF, |<<L,, 
l. = the inner scale of turbulence, 
L. = the outer scale of turbulence, 
figs) - n(zt,)] = difference of refractuom 


index between position) _ @laviciieee 


B. TEMPERATURE STRUCTURE FUNCTION 


The temperature structure function is 
* —_ —. — 
D.(z,-t) = <(T(F) - T(R,)]2>. (eqn 2.4) 


By dimensional analysis, Kolmogorov showed the structure 


function for temperature in the Kolmogorov model as follow: 


Dor) | = Weeds (eqn 2359 
where r = Pe -r, |) and Eee is the temperature structure 
parameter. The logarithm of the structure function, Tapa 
is shown as a function of log | r,- dl in Figure Za 


Equation 2.5 is valid as long as r lies between the 
outer Scale and the inner scale of turbulence. 

For some larger r, r>L, , the temperature fluctuations 
become uncorrelated and the structure function asymptoti- 
cally approaches a limit. For r>om , and assuming homoge- 
neity, <T(r, )>= <T(r,)>, and the temperature structure 


function becomes 


WZ 


Dia (yee <T(r, )>7)- (eqn 2.6) 


<< 


Pauataon 2 .oussnows that the asymptotic limit of D,(r) is 


just twice the variance of the temperature fluctuations. 


fe AMPLITUDE FLUCTUATIONS 


Tatarski [Ref.9] derived the expression of amplitude 
fluctuations based on a perturbation theory. These theoret- 
ical results predict that the amplitude fluctuations at a 
point increaSe without bound with increasing propagation 
distance. For the spherical wave , the log amplitude 


variance 1S 
Ta = 0.124 Negeal en | (eqn 2.7) 


where k is 2R/A and CZ is the index structure parameter. 
The index Structure parameter can be obtained from the 


emperature Structure parameter C? 


2S ale P/T2)«10., (eqn 2.8) 
where P is the pressure in millibars and T is the tempera- 
ture in Kelvin. 

It has been found experimentally [Ref.10] that Equation 
2./ is valid when the fluctuations are small ( Gas0.3 ) and 
that the amplitude fluctuations do not increase beyond 
fa ~ 1.0 with either increasing the path length or 


turbulence strength. 


eS 


G D 
LO t 


SORE 275 


Figure 2.1 Log De as a> bette ci onmos Log [r, - x 


D. INTENSITY SCINTILLATIONS 


We find that the intensity of a beam fluctuates with 
time about its average value <I> in the turbulent medium. 
The intensity fluctuation can be measured with a detector 


recording the intensity of the radiation as a functiongen 


time. 
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Usually we can express the scintillation strength from 
the normalized variance of the intensity fluctuation, as 


for lows , 
eas 
(EGU ee) B-eec(1/1))*2 [)<1/1,>? , (eqn 2.9) 


where I = the intensity at the receiver at time t, 
Le = <1 > 
For the case of weak turbulence, it 1S conventional to 
calculate the variance of the logarithm of the amplitude 
fluctuations rather than of the intensity fluctuations. If 
the log amplitude is normally distributed, the log intensity 


variance can be shown to be 
eel te) <In(i it, 2. ana 


Since the log intensity iS proportional to the square of the 


amplitude, then we have following form; 
(Onion T C,.? )- (eqn 2.11) 


This relation is developed by Tatarski [Ref.9] based on the 


Rytov approximation which gives quiet good agreement for 
we < 0.3. 


E. PATH WEIGHTING FUNCTION 


Tatarski showed that for scintillations the weighting of 
mubpulence along the optical path is not uniform. For the 
case of a reflected spherical wave, the optical scintilla- 
tion weighting function was investigated by Ze'evi in his 
doctoral thesis at NPS. In his work, Dr. Ze'evi assumed 
a plane mirror as_ the target and a weak Kolmogorov 


turbulence model. 
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He derived the weighting function from the log intensasee 
equation of a spherical wave for the folded path using the 
results of Smith [Ref.11]. The log intensity equation for 
the folded path is 


KZ Glee) 
CF C er ei it “k f cz) | Kg¢(K) Sin Re ee (eqn 2. Zam 
where k = wave number ( 2r/A ), 


$(K) = the isotopic form of the Fourier transform 
of the refractive index covariance function, 
ZL = path Voearvom,. 
K = two dimensional spatial wave number, 
of turbulence (= 2i/ L,), 
L = pathlength. 


Ze'evi simplified thiS equation using the Kolmogorov 
Spectrum and the Mellin transform (formula 6.5(15) of 
[Ref.12]) as follow; 


e SL 
2, - ae - (ee 
CiC-0) aan \ nc Naeem a= (eqn 2.13) 
¢ 
where A,= .033(-4)(27°) [(-5/6)cos(5i/12)i?*k”* 
Making the transformation by " X = Z/L ", then we get 
: Jen ae ee (eqn 2.14) 
Gur er) = ALC, \ dxLx(2-K) 


© 


The log intensity for the direct path of length Lugs 


4 OD 


Gu) = ack? | dz| KguK) sin? [EAP 7h hic (edn 2am 


Using the same steps as in the previous work, we get 
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2 eee ves eZ) eu 
Gi= A.Gr(dz(=E] , 
; / 


, ; ah 
where A = .033(-4) [°(-5/6)cos(5@/12)i?k” "= 2 
miking XX = Zj/L ” then 


(ecime. 16) 


é 
ae 


\ 
(ie ALC | AX LX UI-X) (eqn 2.17) 
kK 
Dr. Ze'evi proposed the path weighting functions for the 
Spherical wave by regarding C,; as constance, these are 
Direct path: W, (X) Pere exer) 27* 
Folded path : W.(X) Px x -) 


where X = 2Z/L, 
Z = turbulence location from the detector, 
L = pathlength. 
The two functions are shown in Figure 2.2. AS we can see 


from Figure 2.2, the main contribution for the folded path 


comes from the vicinity of the target mirror. 


ee 
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iit, ESPEREMENTIAL APPROACH 


A. BACKROUND 


Because of several previous attempts to prove Ze'‘evi's 
weighting function for the folded path, most of the equip- 
ment to perform the experiment were available. We made only 
minor modifications to the turbulence chamber that was used 
in the earlier experiments. 

To verify Ze'evi's theory, the turbulence in the turbu- 
lence chamber must be in the Kolmogorov inertial subrange 
region and needs to be as uniform as possible at each of the 
different locations along the beam path. To verify the 
existence of a Kolmogorov structure function, the turbulence 
chamber was first assembled in the basement of Spanagel Hall 
and then tested. After setting the experimental system on 
the roof of Spanagel Hall the turbulence chamber was tested 
again andthe results compared to those measured in the 
basement [Ref.6]. 

During the experiment, we measured Ci? to normalize the 


experimental results at the different locations along the 


beam propagation path. Additionally the optical system had 
to be isolated from the outside disturbances. For this 
reason, we sealed the tunnel with masking tape and the 


optical equipment were also sealed under the tables at the 


end. 


With these considerations in mind, we will discuss the 


entire experimental set up. 
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B. TUNNEL AND TURBULENCE CHAMBER 


We used the same _ tunnel that was constructed ae 
Constantine [Ref.3] and Flenniken [Ref.4]. The tunnel is 6l 
meters long by .6l1 meter wide by .61 meter high and 
consisted of 25 sections 2.4 meters long. The tunnel was 
treated with a water sealant, painted white on the exterior 
and flat black on the interior to eliminate reflection. 


The turbulence chamber was made by modifying a tunnel 


section. It could be placed at any desired location along 
the 61 meter tunnel. Figure 3.1 shows the turbulence 
chamber. A .64 meter high chimney was constructed over the 


opening of the tunnel and a fan was installed in the chimney 
cover to create a positive air flow upward through the beam. 
Above the optical path. a small slit was cut in the side of 
the turbulence chamber to insert the platinum-resistance 
temperature-probes into the region of turbulence. To allow 
airflow in under the heat Source, a wide slit was cut just 
below the heating wire. Both of the holes were baffled to 
prevent unwanted airflow at the turbulence region. in addi - 
tion, two plywood baffles were installed inside the tunnel 
on both sides of the heat Source to reduce the longitudinal 
airflow as much as possible. These baffles had a hole cut 
of sufficient size to allow laser beam propagation without 


obscuring the laser beam. 


C. HEAT SOURCE 


The heat source is shown in Figure 3.2. The heating 
wire was installed on a wooden stand with ceramic insulators 
in ae Serves cof ms Vo Se It was connected to a variac for 
providing the desired voltage. 

The wooden stand was reconstructed so that its width was 
the same as the interior width of the turbulence chamber and 


the height of the. stand was lowered. Additionally we 
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installed a fan rotating at one Hz under the heating wire to 
break the convective laminar flow into turbulent flow. 

The whole process for the turbulence is as follows: The 
exterior air comes into a side slit and is then heated by 
the heating wire. The heated air is broken up into a turbu- 
lent flow by the fan and is drawn upward through the beam 
path by the upper exhaust fan. 


VEPER FAN 





Figure 3.1 Turbulence Chamber. 


Pe TURBULENCE MEASUREMENT EQUIPMENT 


For getting some quantitative measure of the turbulence 
present, the temperature structure constant, Cy? was meas- 
ured at each path position of the turbulence section during 


the experimental run. The C,* was measured using two 
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platinum-resistance temperature-probes separated by 6 centi- 
meter. A schematic of the AC bridge system is shown in 
Fipure 323. 

The resistance of the platinum probes varied in propor- 
tion to the temperature changes in the turbulent’ medium. 
The varying resistance produced the voltage difference in 
the bridge. The fluctuations of the voltage difference in 
the bridge were amplified using a Princeton Applied Research 
model 113 Amplifier then fed to a DATA-6000 waveform 
analyzer. The DATA-6000 then calculated the temperature 


structure parameter. 


ESOP LTIGAL EQULMENE 


We only measured the weighting function for the folded 
path since the direct path scintillation 1S thougnt Seca. 
well explained by the theoms 
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Figure an Schematic of AC Bridge System. 
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The experimental set up for the folded path is shown in 
Figure 3.4. A laser beam from a single mode He-Ne gas laser 
ex 6328A } was used as a source. The laser beam was 
chopped by a mechanical chopper wheel with an open to closed 
rat romwot - eae The chopper wheel rotated at 3000 rpm and 
provided a pulse rate of one KHz. The chopped beam was then 
directed to the target reflector by a beam spliter. ; 

We employed a corner cube, flat mirror and a cat's eye 
as the target reflector and concentrated on the comparison 
of the scintillation results from the different reflectors. 

The optical scintillation was detected by an RCA #30827 
Silicon Avalanche Photodiode. When operated at 250 volts, 
the detector provided extremely high signal to noise ratio 
without the need for cooling. We did not uSe any spectral 
filters as all of our experimental work was done inside the 
tunnel at night. . 

The output from the detector was amplified by a PAR 
model 113 amplifier. The high and low frequency roll-offs 
were set to allow the signal to be amplified without 


creating any distortions of the signal. 


F. DATA PROCESSING EQUIPMENT 


The signals amplified by the PAR #113 were passed to the 
demodulator. This device sampled an incoming signal and 
held its instantaneous maximum signal level for the desired 
duFarron. The background was then sampled and the differ- 
ence of these two levels was taken. This difference was 
held constant until the next difference was measured. The 
reconstructed signal was then passed to a log convertor. 

The log-convertor produced an output voltage that was 


proportional to the lop of the tnputevotezasee 
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G. DATA REDUCTION SYSTEM 


ite wcimeput sar ethe lop ‘convertergmerthe clock pulse of the 
chopper wheel, and the amplified temperature fluctuation 
Signal were transfered to the DATA-6000 waveform analyzer. 
We used a Model-610 plug-in module which has sampling rates 
up to 100,000 per’ second. AD SHEA SS 254A" catculator was 
programmed to calculate the standard deviation of the ampli- 
fied temperature fluctuations and of the log intensities. 

We made five separate runs at each turbulence chamber 
location and each run consisted of 500 samples. The temper- 
were Structure constant, C,*, and the log intensity vari- 
ance were then calculated from these data by the HP 9825A 


ealculator. 
The HP 9825A calculator functioned as the controller and 


input-output device to the DATA-6000, He 7coce plotter, 
and HP 9871A printer. 
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Figure 3.4 Experimental Set Up for the Folded Path. 
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Vee Ch Ree NT AG WORK 


A. GOALS 


There were three main goals in our experiment. We 
continued essentially the same experiment as the work of 
Constantine [Ref.3] and Henry [Ref.5] with the following 
Seals in mind. 

First, to create turbulence in a controlled region that 
approximated the Kolmogorov spectrum. 

Second, to prove or disprove Ze'evi's weighting function 
for the folded path experimentally. 

ard , to compare the results uSing three different 
finds of reflectors ( corner cube, flat mirror, cat's eye ). 

The corner cube had been chosen by Henry because of it's 
complicated reflection pattern. Additionally, we employed 


the cat's eye because of it's low reflection. 


eee PROCEDURE 


In order to achieve Kolmogorov turbulence in the 
chamber, many runs . were conducted in the basement of 
Spanagel Hall. After assembling the whole experimental 


system on the roof of Spanagel Hall, we retested the turbu- 
lence chamber and obtained the same results that we obtained 
in the basement to verify the region of the Kolmogorov iner- 
tial subrange. 

All optical equipment were taped to the concrete floor 
and the reflector was on an optical bench which was anchored 
with lead bricks for stability. The tunnel sections were 


also anchored with lead bricks for wind stability. 


ZT 


We performed the experiments after dark in order to 
reduce the turbulence in the tunnel caused by the _ sun 
heating the tunnel. This also reduced the noise caused by 
light leaks in the tunnel. The tunnel sections were taped 
to prevent wind disturbances. 

The quiet tunnel scintillation was measured for a refer- 
ence before we turned on the heater. The scintillation and 
the variance of temperature fluctuations were measured 
Simultaneously by the DATA-6000 waveform analyzer. The 
average temperature inside the tunnel was measSureed by a 
copper constantan thermocouple using an ice water reference. 
The turbulence was then relocated along the path and 
retaped. This procedure was repeated for each data point. 

The DATA-6000 calculated the log intensity variance 
which was then recorded on magnetic tape by the HP 9825A 
Caleu lator. 

The log intensity variance of the quiet tunnel was 
Subtracted from the log intensity variance at each point. 
This difference of the variances was then normalized by a 
Factor of | Gee-G.*> Gree ae The T and C,* are the 
temperature and temperature structure constant at the point 
respectively, < > is the ensemble average of the variable 
over the path. These normalized variances were plotted 


versus path position. 


Corer ER IME 


In this section, we discuss briefly the experimental 
measurements discussed in detail by Mr. Kim [Refs7am 
present the experimental details in his paper. We employed 
three different kinds of reflectors, (corner cube, flat 
MacteisOras and cat's eye), performed experiments for each 


reflector, and compared the results each other. We first 
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tried to measure the scintillation for the three reflectors 


during the same runs. Dye Lacincmeme trek bectors on an 
optical bench, we could change a reflector quickly. During 
the experiment, the alignment of the flat mirror was 


destroyed when we changed to another reflector even though 
the change was made carefully. Additionally we observed an 
interference pattern on the detector face. So we attached a 
small black tape to intercept the secondary beam off the 
back face of the splitter. Finally, we decided to make 
measurements on the reflectors separately. Because of the 
problem of interference from multiple reflections, we made a 
Screen in the detector's back board to observe the reflected 
beam structure at the detector. Since intensity saturation 
was also observed we decided to employ a convex lens between 
the source and the chopper. We also observed GZ > .3, so we 


reduced the power to the heater in the turbulence chamber. 


Because of outdoor work, the experiments had to’ be taken 
during a calm night. We choose the turbulence location at 
roe ley 74 eee Ss . 06, 30750,°-3/7.70, 42.58, 47.46, 52.34, 
and 57.22 meters along the path _ position. The data 


collected involved 5 sets, of 2500 samples each for each 
turbulence location. The data was then processed _ and 


plotted as discussed in the previous section. 


lL. Corner Cube 
We took data on ll Oct, ZO Geren uaa 17 sOct . The 
lower power was used for the heater. We retaped the whole 
tunnel before making each experimental run. After the first 


run we added the baffles at both sides of the heat source 
and at the slits of the turbulence chamber to limit unwanted 
drafts through the turbulence chamber. Figure 4.1, Figure 
4.2, Figure 4.3 show the plot of the results. The scintilla- 
tions had an unexpected increase near the reflector end of 


the tunnel. The data are listed in Appendix A. 
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2 Flat Mirror 


The experiments for the flat mirror were taken on 20 
Oct, 22 Oct, and 24 Oct at night. The weather was clear with 
a light breeze. 

We rechecked all of the equipment andthe tunnel 
section but couldn't find any problems. The scintillation 
still increased near the vicinity of the reflector but more 
strongly than the corner cube as shown in Figure 4.4, Figure 


4.5, Figure 4.6, and the data in Appendix B. 


3. Gat s Eye 


It took many trials to align the cat's eye because 
the beam returned from the Cat's eye appeared to have 
interference-fringe bands on the order of 5 mm wide at the 
detector plane, and the intensity of the returned beam was 
low. Therefore we removed the convex lens from the front of 
the laser to increase the back scattered laser energy and 
made measurements with the distorted beam Seen During 
the first run for the cat's eye, we realized the detector 
and the source were tilted when we moved the turbulence 
chamber. Therefore we measured scintillation from the 15.74 
meter point for the first run. Since the scintillation data 
for the first wo experiments were too low, we increased the 
power to the heater. 

As we can see in Figure 4.7, Figure 4.8, Figure 4.9, 
the scintillation gradually dropped from the 30.38 meter 
point to the end position. The detailed results are in 


Appendix C. 
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V. DISCUSSION AND OBSERVATIONS 


Loaceehas Chapter. we consider the size of the 
Fresnel-zone-eddies and then measured the image point of the 
eddy lens combinations from the detector due to the turbu- 
lent eddy location along the path length. 

Finally we discuss the experimental results in a view of 


the above discussed points. 


A. GEOMETRICAL MODEL OF TURBULENCE 


The physical mechanism by which turbulence induces beam 
break-up is random focusing of light by refractive inhomo- 
genelities in the path of the optical wave. The geometrical 
expressions were first developed by Tatarski [Ref.9]. He 
described the eddies in the atmosphere as spherical lenses 
whose refractive index is determined by the fluctuation dn, 
from their Surroundings. | 

Since dn is much less than unity, we can neglect the 
reflection at the surface of the random medium. Therefore, 
we did not consider the effects of absorption or scattering 


by either aerosol or molecules in this paper. 


In tne paraxial approximation, it is well known that an 
eddy can act as a lens with a focal length f = r/dn (where 
r is the radius of eddy). For dn>O, the lens is convergent, 
and for dn<O, the lens is divergent. In the atmosphere, the 
fluctuation dn is much less than Mo, which implies that 
£200 3K: 


Let's consider the turbulent eddy of radius r in the 
atmosphere as Shown in Figure 5.1. The spherical wave of 
a wavelength NX encounters a_ turbulent eddy at position A 


between source S and the observer's plane P. The eddy that 
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Figure 5.1 The Geometry of the Simple Turbulent Eddy. 


is the most effective in producing Scintillations at the 
observer's plane is the one that satisfies the condition 
SBP - SP = A/2, so that the two rays interfere destruc- 
tively. Working out the geometry of Figure 5.1, Clifford et 
al [Ref.13] showed that this Fresnel-zone-size eddy has a 


radius 
meme ZA(l =e Z / i) (eqn 5.1) 


where Z is the distance of the eddy lens from the source and 
L is the path length. 

Let us consider the effects of the various size eddies 
more the scintillation. Eddies smaller than Fresnel zone 
contribute less because of the weak refractive fluctuations. 


Larger eddies also will not’ produce strong Sombie lat POnwat 
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the observer plane because of their long focal lengths. 
Even the smallest details in the scintillation pattern 
arrive at P unaffected, that 18, unsmeared by the large- 
scale distortions. 

The sizes of the Fresnel-zone-eddy along the path length 
of 61 meters for the direct and for the folded path were 
plotted in Figure 5.2 for a wavelength ot ysa5Z5o The 
calculated values are listed in Appendix D. 

If we assume that the turbulence has a Kolmogorov spec- 
trum, it is clear that the variance of the refractive-index 
fluctuations attributable to irregularities of optimum size 
varies systematically along the _ path. Consequently the 
scintillation weighting function is a function of loeataiae 
along the path 

From Figure 5.2 it is clear that this weighting funemeom 
must reach a maximum at the midpoint of the path and for the 
direct path it has to drop symmetrically to Zero at jeme 
ends. 

Furthermore, the Scintillation weighting func vom 
assumes that the turbulence has a Kolmogorov spectrum. Lt 
is possible for the size of the most effective eddies to be 


smaller (mm) than the Kolmogorov inner scale. 


B. THE FOCAL POINTS OF SINGLE EDDY FOR Dib Orr ny ee ia 


Suppose that the turbulent eddies are two positive 
lenses of the same focal length (f) and at the distance X 
from each end. This can be considered as a combination with 
two turbulent eddy lenses on the direct path. The schematic 
geometry is shown in Figure 5.3 This model has’ been over 
Simplified by use of an eddy lens model. (Actually it would 
be more rigorous to use a set of orthonormal functions to 


describe the turbulent medium. ) 
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Figure 5.2 The Radius of the Eddy of Fresnel-Zone 
—EewtOEeDmnect Pati (roa) and for 
Folded Path Cc) along a 61m Path. 


From the Gaussian lens equation the focal points of the 


eddy lens combination can be obtained. 


— 2fUL-X)- £X/(K- f) (eqn 5.2) 
: L(L-X)-f- TA/(x-6) 
where Ss, = ilmage point of the eddy lens combination, 

= focal length of the eddy, 


f 
L = path length ( source to reflector ), 
X 


= distance from the source. 
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From equation 5.2, the image distances from the receiver 


plane R is 


als Si aN (eqn 5 aaa 
Taking X' = X/L, £° = £/L then equation Se3ebeceomes 
= ss;2d 97 « 
: 2 UR KI- F KS OK FY , 
(= ees ee (eqn 5.4) 


21- xD = Fe FKIR FD 


The numerical values along with the various focal length are 
listed in Appendix E. 

As we can see from the Appendix E, when f < L, the beam 
focuses at the receiver plane when the turbulent eddy is 
located at xX = f. AG: <f7= ee the beam focuses’ at the 
receiver plane when the turbulent eddy lens is at  “ehe 
mirror, and with a longer focal lengths than path length L, 
the beam focuses beyond the receiver plane, if at alm 
Therefore, in order to focus at the receiver plane the 
focal length (f) of the eddies must be less than path 
lene tin. (bE) 


C. BEAM SPREADING 


If the laser beam propagates through the random medium, 
the phase fluctuations in the wavefront will produce a meas- 
urable spreading of the beam. The analysis of this effect is 
SUEpriSine ] y “di Prmeniie. though an approximate geometrical 
solution has been presented by Yura [Ref.14]. Lawrence 
[Ref.10] explained its approximate nature from the following 
simple argument. Referring to Figure 5.4, consider a laser 
source § emitting a well collimated laser of a spot diameter 
don the observer plane R. If the parcel of turbulence were 


inserted in the beam at point A, the beam spreading would 
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Fogure 5.3 The Schematic Geometry of the Two Lens Combination. 


arise with angle © and make a spot of diameter L at observer 
plane. et. instead, the same parcel of turbulence is 
located at point B, the beam spreading caused by it will 
still have the same angle 6, but the spot on the receiving 
screen will be only dp in diameter. Thus, if the inherent 
angular divergence of the laser beam is much’ smaller than 
the atmospheric effect, the weighting Suneciaom that 
describes the effectiveness for beam spreading of turbulence 
at various positions along the path is'7 linear. It ais a 
maximum at the source and drops to zero at the observer 


plane. 
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beam Source 





Figure 5.4 The Geometric Figure of Beam Spreading 
at the Observer Plane. 


D. DISCUSSION AND CONCLUSIONS 


We have mentioned that the Fresnel-zone-size turbulent 
eddies along the path can mainly contribute to the scintil- 
lation at the detector. That 2s the size of the 
Fresnel-zone-eddy is determined by the path length and by 
the path position as we saw in equation 5.1 

As we can see from the Appendix D, the Fresnel-zone-size 
eddies have the diameter of 1.4 mm up to 6.6mm _= for the 
direct path and of 1.6 mm up to 8.8 mm_ for the folded path 
of the 61 meter tunnel. 

| From Mr.Lee's works we realized that the inner scale of 


turbulence in our tunnel is 7.5 mm and the outer scale is 
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67.5 mm. These eddies have much larger diameter than the 
Fresnel-zone eddies for the 61 meter tunnel. eras clear 
Ehat these eddies do not produce strong scintillation in the 
short path position as we discussed in the section 5.1. 

In our experimental results, most of the measured scin- 
tillation from 1.1m to the midpoint of the tunnel showed 
lower values of scintillation than the theoretical values as 
we can see from the figures in chapter 4. 

The strong scintillation results near the reflector can 
be explained as follows: 

First, the size of the Fresnel-zone at the 30 m position 
is around 7.5 mm. This is of the order of the inner scale of 
our turbulence. Therefore, from the 30 m path position, we 
can consider that Fresnel-size-zone eddies exist along the 
math position and they contributed to the scintillation at 
the detector. Second, we must also consider the intensity of 
the reflected beam. As we discussed in the previous section, 
the beam spreading decreases when the parcel of turbulence 
is located near the reflector. When the reflector is 
smaller than the spot size, some portion of the total inten- 
Sity should be lost. Since the intensity profile at the 
reflector has bright and dark regions on the order of one 
Fresnel-zone in diameter and these regions are moving with 
the turbulent flow, a small reflector will not be fully 
1lliminated. This will induce additional scintillations in 
the reflected beam. 

Having discussed a simple geometric model of turbulence 
and analysis of our experimental results, we conclude that 
we are not able to verify the Ze'evi's weighting function 
for the folded path with a short path length because of the 
absence of Fresnel-zone-size eddies. 

We recomand that the tunnel should be longer to verify 
the Ze'evi's theory and the reflector must be larger than 


the size of beam spot. 
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APPENDIX A 
CORNER CUBE 


DATA OF THE FIRST EXPERIMENT, El) OCT aia 


TURBULENCE NORMALIZED 
LOCATION Ct (k-7 n- a ee LOG- INTENSITY 
QUITE 1227 + .0039 
8.42 25.41 .1488 + .0013 1.106 + .058 
15.74 o5eo0 .1820 + .0015 2.176 + .052 
23.06 31.91 .2011 + .0060 3.205 + .102 
30258 3877. 2434 + .0034 4.558 + .143 
WAS TAS. 81.73 .2721 + .0019 2.863 + .085 
42.58 Be 22 .2901 + .0071 7.054 + .135 
47.46 61.03 3363 + .0062 5.925 + .152 
52.34 28.54 .2825 + .0019 8.797 + .409 
57.22 9.77 2296 + .0073 15.160 + .124 
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pane OFS THE SECONDMEXPERIMENT, 12. OCT. 1984 


TURBULENCE NORMALIZED 
LOCATION Ct?(K?2/m? ) iy & LOG- INTENSITY 
QUITE .0477 + .0002 
1 .uG 20 ia 0629 + .0021 355 + .001 
8.42 24.25 .1245 + .0219 DEO, 1911 
15.74 52.52 .1841 + .0025 2.624 + .130 
23.06 58.54 2278 + .0074 3.663 + .206 
30.38 67.69 2954 + .0145 5.446 + .217 
37.70 48.99 27 + . 0093 6.503 + .227 
42.58 52.19 .2954 + .0005 7.165 + .159 
47.46 53.24 .3014 + .0083 7.343 + .233 
52.34 44.03 3005 + .0100 8.993 + .374 
57.22 15.97 2576 + .0146 17.770r + .329 
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DATA OF THE THIRD EXPERIMENT, l723CCl ee 


TURBULENCE NORMALIZED 
LOCATION Ct?(K?/m? ) 2s LOG - INTENSITY 
QUITE .0577 + .0050 
1.10 9.66 0607 + 0022 .088 + .009 
8.42 12.51 1267 + .00@8 2.535 + .115 
15.74 57.86 .1840 + .0057 1.319 + .162 
23.06 La G7 .2390 + .0104 3.200 + .145 
30.38 19.61 .2668 + .0026 8.385 + .16l 
37.70 270 .2763 + .0106 6.569 + .338 
42.58 10.90 227 Gn OO 10.603 + .056 
47.46 25 37 .3290 + .0070 9.810 + .281 
52.34 10.55 DRE. = Oome 16.100 + .116 
S722 28.08 3141 + .0158 7.897 + .416 
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APPENDIX B 
FLAT MIRROR 


DATA OF THE FOURTH EXPERIMENT, 20. OCT. 1984 


TURBULENCE NORMALIZED 
LOCATION Ct2(K?/m?_ ) 1G LOG- INTENSITY 
QUITE .1887 + .0010 
6 roy .1907 + .0005 238 + .0O11 
gi 42 66.35 .2141 + .0005 .566 + .025 
oe 40.84 2529 + .0001 2.515 + .181 
23.06 35 9a" .2914 + .0049 4.888 + .268 
30.38 94.59 3770 + . om 4.051 + .416 
S70 42.96 .4376 + .0011 12.980 + .489 
42.58 17 a7 4292 + .0156 30 = 2c Ome eal 
47.46 20.89 .4877 + .0187 34.780 + .659 
2 34 levee 5246 + .0088 49.070 + .615 
57 oe 7.60 .4830 + .0527 91.460+ 1.270 
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DATA OF THE FIFTH EXPERIMENT, 225 5@C Ree 


TURBULENCE NORMALIZED 
LOCATION Ct?(K?/m? ) <i LOG- INTENSITY 
OUT E .1868 + .0013 
ecu 9.17 1921 + .0008 .698 + .018 
8.42 3 oe .2218 + .0010 1.199 + .068 
15.74 38.84 2649 + .0099 2.842 + .189 
23.06 50.16 - .3093 + .0005 3.715 + .166 
30.38 LAB .3784 + .0016 6.932 + .237 
37.70 42.03 .4382 + .0049 11.400 + .162 
42.58 330s .4930 + .0138 19.280 + .926 
47.46 24.98 4445 + .0023 20.180 + .513 
52.34 14.34 4651 + .0101 38.690 + .800 
Sn oe 11.19 4540 + .0070 46.530 + .753 


Do 


ao 


DATA OF THE SIXTH EXPERIMENT, 24. OCT. 1984 
TURBULENCE NORMALIZED 
ROC ROUEN Cl (K-/m- ) < J > LOG- INTENSITY 
QUITE wey > + 0005 
ite 14.99 .1884 + .0008 1.448 + .017 
8.42 19.98 .2165 + .0047 Peceee 5072 
15.74 Oe .2487 + .0060 3.636 + .128 
23.06 28.47 .2984 + .0088 6.146 + .161 
30.38 23 5c7 eee 0272 10.602 + .496 
S77 0 46.95 eames sOil0 8.224 + .602 
42.58 28.98 .4178 + .0062 14.639 + .464 
47.46 5 Onmay .4698 + .0065 18.183 + .889 
52.34 DAL. OD .4788 + .0080 26.879 + .669 
Se 39.93 4951 + .0184 IS) D7 i ean 
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Gets SHE 


DATA OF THE SEVENTH EXPERIMENT (ZG yee 


TURBULENCE NORMALIZED 

VOCAT HON eet (K-71 am enitue LOG- INTENSITY 
QUITE .0779 + .0003 
15.74 39.36 .1510 + .0049 1.003 + .044 
23.06 31.26 .1701 + .0083 1.709 + .125 
30.38 30ne2 .2120 + .0016 2.231 + .086 
37270 24 Sie 2087 = Oem 3.206 + .227 
42.58 15.48 l701 = . 00Rm 3.351 + .089 
47.46 11.50 .1705 + .0078 4.454 + .140 
52.34 10.94 .1336 + .0049 2.394 + .065 
57.22 8.79 .1339 + .0021 3.001 + .050 
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DAP wOPwiHr ELGHI BN EXPERIMENT, &. NOV, 1984 


TURBULENCE NORMALIZED 
LOCATION Ct?(K?/m? ) 21S LOG- INTENSITY 
QUITE .1564 + .0029 
3.55 30.25 .1703 + .0042 .602 + .086 
8.42 2038 .1810 + .0013 1.092 + .027 
15.74 55.25 .2118 + .0005 Haas 037 
DeneG 44.28 — .2431 + .0005 3.078 + .190 
30.38 53.09 miese + 0010 3.286 + .085 
7710 55.78 .2844 + .0029 3.945 + .107 
42.58 47.25 .2714 + .0091 4.152 + .206 
47.46 40.44 .2724 + .0010 4.882 + .067 
52.34 24.84 .2543 + .0003 6.498 + .254 
57.22 12.96 .1956 + .0083 4.248 + .113 
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DATA OF THE NINTH EXPERIMENT, Ll NOV eee 


TURBULENCE NORMALIZED 
LOCATION ) Ce*(K?2/m? ) - LOG- INTENSITY 
QUITE OO Se COS 
35D 2 38402 HL250 0) O02) 471 2 ie 
8.42 122s 1708) 0055 1.999 + J1Se 
15.74 sie) sue 2360, = -0208 3.975 2 ame 
23506 101 07am .2886 + .0010 5.912 fee 
30.38 OS OS 23028 + &CGza 7.974 2 “Zep 
Be710 Somes 2993) 2 00cm 14.425 + aye 
42.58 53.43 1250167 => = OnelO 11.390 + .204 
47.46 34.64 2501, + J0148 16.122 4) eae 
52.34 64.26 2072 + 2 Oden 5.741 = eige 
S72 6 Loy “i803 + Ope 4.385) ee 
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APPENDIX D 
FRESNEL-ZONE-SIZE EDDIE'S RADIUS FOR A 61 METER PATH 


TURBULENCE EDDY RADIUS (millimeter) 
LOCATION DeREGl PATH BOLDED fh aE 
SOURCE 
Ihe Ue) eC a7, 7o 30 
8.42 2.143 ted 
oe 2./18 2.945 
2 100 beO 3.440 
SO ARSIS: D. LOG 3.600 
By. 7 0 or O18 4.060 
42.58 ZO 53 4.188 
47.46 ZS 2 4.284 
O12). 34 Zor 4.349 
Dy a2 2 1.498 4.385 
60.00 Og ae 2 
oll OO .00 oS 9S 
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APPENDIX E 
POSITION AT FOCAL POINTS MEASURED FROM THE DETECTOR 


EDDY LENS 
POSITION FOCAL POINTS (aL) 

FROM DETECTOR f=.1L f=.20, £=.31 £ = 94 
.00 O.11 0.22 0.35 0.50 0.60 0. 
be 0.00°— 0: 12uOnes 0.40 0.58 0. 
aie -0.09 0.00 0.15 0.30 0.47 0. 
ci =-0, 19° =0. 030 OeOG 0.17 0.35 OF 
one -0.29 -0.13 -0.40 0.00 0.19 0. 
.5L -0539 -0.21) 052200 ) oan 0.00 0. 
Rol 0.48 =0,27 Weet20 Sie oe 0.19 On 
Pile -0.57. -0.17 -0.80 -0.52 =0.35) eicE 
ue -0.65 -1.20 -0.74 -=0.60 9 =0047.-GF 
OL -1.01 -0.86 -0.77 -=0.69 =0259) oie 

Reflector -0.95 -0.89 -0.75 -0.67 -0.57 -0O. 
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EDDY LENS 


POSITION FOCAL POINTS (4L) 

FROM DETECTOR f=.7L f=.8L f£=.9L f= L  f=1.1L £f=1.2L 
00 I ORIN. 64 2.00 2.44 3.00 
a CO RenRIee? 7 1 GO 1.99 2.49 3.13 
m2 Oo me 201s 55 1.98 2.53 3.25 
at OC ORO Sm aS Le Sil 2.50 a0 
AL On) SO: IX 1.75 Dave 3.20 
5L 0.44 0.72 1.06 1.50 2.08 2.89 
fan 0.21 0.47 0.78 eet ae 2.40 
a7 Of00 7 0.22 0.48 0.82 1.25 1.83 
8L TOR ONO CMO 2? 0.49 0.84 eo 
core FO seme 118 Oe OO 0.22 0.50 0.85 

Reflector -0.46 -0.33 -0.18 0.00 2? 0.50 
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EDDY LENS 
POSITION 


FROM DETECTOR £=1.3L £=1.4L 
4. 
16 
72. 
we 
44 
09 
16 
97 
86 
98 
Boe 


00 
els 
<216 
3 
~4L 
ois 
.6L 
ye 
.8L 
ae 
Reflector 


oe 


wi 


3.98 
“os 
4.45 
4.44 
4. 
3 
2 
1 
1 
0 


10 


45 
207 
aoe 
3 
206 
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FOCAL POINTS (AL) 
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